Search for B ->■ h<*W with the full Belle T(4S) data sample 
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We report a search for the rare decays B — > h,(*'vp, where h^ stands for K + , K s , K* + , K*°, 7r + , 
7t°, p + , p° and <f). The results are obtained from a 711 fb" 1 data sample that contains 772 x W 6 BB 
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search for signal candidates by fully reconstructing a hadronic decay of the accompanying B meson 
and requiring a sing le /i w meson left on the signal side. No significant signal is observed and we set 
upper limits on the branching fractions at 90% confidence level. The measurements of B + — > K* + vv, 
B + — > ty + uu, B° — > tPvv and B° — > p°vv provide the world's currently most restrictive limits. 
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The decays B — > K^vv proceed through the 
flavor-changing neutral-current process b — > svv, which 
is sensitive to physics beyond the standard model 
(SM) [TJ The dominant SM diagrams are shown 
in Fig. [TJ The SM branching fractions are estimated 
to be (6.8 ± 2.0) x 1CT 6 for B+ ->■ K* + uu [2] and 
(4.4 ± 1.5) x 1(T 6 for B+ -> K+vv decays [T]. The 
decays B — > (jr, p)vv proceed similarly through b — > dvv. 
Compared to b — > svv transitions, the branching 
fractions are further suppressed by a factor | Vtd/Vts | 2 - 
The decay B° — > 4>vv proceeds through a yet unob- 
served penguin annihilation process, with the expected 
branching fraction thus much lower. The advantage of 
vv rather than £ + i~ in the final state is the absence of 
long-distance electromagnetic interactions. In the ratio 
of the individual branching fractions for B — > Kvv and 
B — > Kl + t~ , the form factor normalization cancels out, 
leading to a factor of three smaller theoretical error 



compared to the vv mode alone [2]. Measurements of 
the B — > Kvv and B — > K£ + £~ branching fractions 
might reveal moderate deviations from SM expectations 
due to New physics such as SUSY particles, a possible 
fourth generation and a non-standard Z-coupling which 
would contribute to the penguin loop or box diagram 
and affect the branching fractions [3]. 

Experimental measurements [4] of the b — > s£ + £~ 
transitions with two charged leptons are in good 
agreement with SM calculations [2J. The challenging 
search for decays with two final-state neutrinos was 
previously carried out by the CLEO, BaBar and Belle 
collaborations [5-8 . No signal was observed, and the 
experimental upper limit for the B + — > K + vv decay is 
a factor of three above the SM prediction; for the other 
branching fractions, the limits are an order of magnitude 
above the predictions. 
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This measurement of B — > h^fv, where ftW stands 
for K+, K*+, K*°, tt+, tt , p+, p° and (j> [§], is based 
on the full Belle data sample recorded at the T(AS) 
resonance that contains 772 x 1Q 6 BB pairs. The main 
improvements compared to the previous analysis [5] 
consist of the use of a new probabilistic full reconstruc- 
tion, a further optimized background suppression and 
an improved signal extraction procedure. 



b 




(a) Penguin diagram (b) Box diagram 

FIG. 1: The quark-level diagrams for the b — > svv transition 
in the standard model. 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux return located outside of 
the coil is instrumented to detect K\ mesons and to 
identify muons (KLM). The detector is described in 
detail elsewhere [ID] . Two inner detector configurations 
were used. A 2.0 cm radius beam pipe and a 3-layer 
silicon vertex detector was used for the first sample of 
152 x 10 6 -BB pairs, while a 1.5 cm radius beam pipe, 
a 4-layer silicon detector and a small-cell inner drift 
chamber were used to record the remaining 620 x 10 6 BB 
pairs [TT] . The data set recorded with the second config- 
uration of the SVD was reprocessed with respect to [TO] 
using new track finding algorithms, which improved the 
track reconstruction efficiency. A GEANT3-based [12] 
Monte Carlo (MC) simulation of the Belle detector is 
used to optimize the event selection and to estimate the 
signal efficiency. 

We identify signal candidates by fully reconstructing 
the accompanying B meson (Stag) an d requiring one 
single meson on the signal side. The B tag candidates 
are reconstructed in hadronic decay channels using a 
neural network-based hierarchical full reconstruction 
method [13] . which provides, for a given purity, roughly 
twice as many B tag candidates compared to the full 
reconstruction method used in the previous analysis [BJ. 
The reconstruction is done in four stages; at each stage, 
the signal probabilities are calculated. In the first stage 
charged tracks, photons and K® and ir° mesons are 
reconstructed. In the following step, two to five of these 



particles are combined in different modes to form -D^, 
D° and J/ip candidates. Some of the most important 
variables used in the neural network training are the 
product of the neural network output for the children, 
the invariant mass of child pairs and the angle between 
them, the angle between the momentum of the D meson 
and the vector between the D decay vertex and the inter- 
action point, and the significance of his vector's length. 
In the third stage, the particles from the prior stages 

are combined to form the DJ^X and D*° mesons. In the 

(») 

final stage, the B^- and B° candidates are reconstructed 
in one of 1104 exclusive hadronic decay channels. 
Here, variables with good discrimination power are the 
product of the neural network outputs of the children, 
the mass of the D meson, the mass difference of the D 
and D* candidates, the angle between the B meson and 
the thrust axis, and angles between the children. We 
use the output variable o tag of the full reconstruction, 
which ranges from zero for background events to unity 
if a clear -B tag is obtained from the network, and require 
otag > 0.02. This cut was found to give the best 
expected branching fraction limit for all channels. We 
select the B tag candidates using the energy difference 
AE = Eb — -Ebcam and the beam-energy constrained 
mass M hc = v^bcam = Pb > where E beam is the beam 
energy and Eb and pb are the reconstructed energy and 
momentum of -Btag candidate in the T(AS) center-of- 
mass (cm.) frame. We require -Btag candidates to satisfy 
the requirements M bc > 5.27 GeV/c 2 and -0.08 GeV 
< AE < 0.06 GeV. If there are multiple B tag candidates 
in an event, the candidate with the highest o tag is chosen. 

The particles in the event not associated with the 
Btag meson are used to reconstruct a B sig — > h^vv 
candidate. Prompt charged tracks are required to have a 
maximum distance to the interaction point (IP) of 5 cm 
in the beam direction (z), of 2 cm in the transverse plane 
(r — (f>) and a minimum momentum of 0.1 GeV/c in the 
transverse plane. K (tt*) candidates are reconstructed 
from charged tracks having a kaon likelihood greater 
than 0.6 (less than 0.4). The kaon likelihood is defined 
by TZk = £-k/(£k + £-k), where Lk{£-k) denotes a 
combined likelihood measurement from the ACC, the 
TOF, and dE/dx from the CDC for the (tt*) tracks. 
It is a function of the momentum and the polar angle 
of the tracks in the laboratory frame. The kaon (pion) 
identification efficiency is 88%-93% (86%-93%) with a 
pion (kaon) misidentification probability of 10%-12% 
(8%-ll%). We use pairs of oppositely charged tracks 
to reconstruct Kg decays, with an invariant mass that 
is within ±15 MeV/c 2 of the nominal K® meson mass 
(corresponding to a width of 5.8cr). The 7r + 7r _ vertex 
is required to be displaced from the IP by a minimum 
distance of 0.22 cm. The direction of the pion pair 
momentum in the transverse plane must agree with the 
direction defined by the IP and the vertex displacement 
within 30 mrad. For ir° — >• 77, a minimum photon 
energy of 50 MeV is required and the 77 invariant mass 
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must be within ±16 MeV/c 2 of the nominal 7r° mass 
(4.1a). 



The decays B+ ig -» 
B®-,„ — > K%vv and £?°- — > 7r u z^ are reconstructed 

falg O big 

from single K + , 7r + , Jfg and 7r° candidates, respectively. 
The £ s ° ig ->■ A*<W candidates are reconstructed from 
a charged pion and an oppositely charged kaon, while 
Bgtg — > A"* + ^P decays are reconstructed from a ATj 
candidate and a charged pion, or a 7r° candidate and a 
charged kaon. The reconstructed mass of the K*°(K* + ) 
candidate is required to be within a ±75 MeV/c 2 window 
around the nominal K*°(K* + ) mass. Furthermore, pairs 
of charged pions with opposite charge are used to form 
B® ig p Q vv candidates, where the tt + tt~ invariant mass 
must be within ±150 MeV/c 2 of the nominal p° mass. 
For figjg — > p^vv, a charged pion and a ir° candidate 
within a ±150 MeV/c 2 mass window around the nominal 
p + mass are used. A K + K~ pair with a reconstructed 
mass within ±10 MeV/c 2 of the nominal <j> mass is used 
to reconstruct <f> candidates. After identifying the Stag 
candidate and reconstructing the light meson, we require 
that no additional charged tracks nor 7r° candidates 
remain in the event. 



The dominant backgrounds are from e + e~ — > qq 
(q = u,d,s,c) continuum events and BB decays with 
a b — > c transition. During the full reconstruction, a 
continuum suppression algorithm based on modified 
Fox-Wolfram moments [13] is applied. To further 
suppress the continuum background, we use the cosine 
of the angle between the momentum of the and 
the thrust axis with the sign convention to the side 
of momentum flow of the rest of the charged tracks, 
evaluated in the T(45) rest frame. This cosine is close to 
— 1 or 1 for continuum events but uniformly distributed 
for spherical BB events. We require the cosine to lie 
between —0.8 and 0.7. In this way, the continuum 
background component is nearly completely removed 
from the signal region, which leads to a better signal 
sensitivity compared to the previous analysis [5J . 

We introduce a lower bound of 1.6 GeV/c on the 
momentum of the M*-* candidate in the B s [ g rest frame 
to suppress the background from b — > c transitions. An 
upper bound of 2.5 GeV/c rejects the contributions 
from radiative two-body modes such as B — > K*j. The 
momentum requirement is removed for <j) candidates 
due to the lack of theoretical calculations for B — s- 4> 
form factors. To suppress backgrounds with undetected 
particles produced along the beam pipe, we require the 
cosine of the angle between the missing momentum in 
the laboratory frame and the beam to lie between —0.86 
and 0.95. Contributions from rare B decays involving 
b — > u, b — > s, or b — > d processes are found to be 
small according to MC studies. The only exception 
is the B s ; g — > 4>uu decay, where rare decays represent 



the majority of the remaining background events. The 
B + — > t + v t decay with the it + vv and p + vD final states 
contributes only 3% and 2% of the total background in 
these channels, respectively. 

The efficiency of the full reconstruction differs be- 
tween data and MC simulation. The correction ratio, 
depending on the Stag decay mode and obtained from a 
study using b — > c semileptonic decays on the signal side, 
lies between 0.7 and 0.8 and is applied to all correctly 
reconstructed -Btag candidates in the MC simulation. 

The most powerful variable to identify the signal 
decays is the residual energy in the ECL, -Eecl, which 
is the sum of the energies of ECL clusters that are 
not associated with the Stag daughters nor with the 
signal-side fv-*' candidate. To suppress contributions 
from noise in the calorimeter, minimum energy thresh- 
olds are required: 50 MeV for the barrel, 100 MeV 
for the forward endcap and 150 MeV for the backward 
endcap region. In a properly reconstructed signal event, 
no activity should appear in the calorimeter, so signal 
events peak at low -Eecl values. 

The signal yield is extracted from an extended binned 
maximum likelihood fit to the -Eecl distribution in the 
range from to 1.2 GeV. The likelihood is 



C 



N 



N 



(i) 



»=i j 



where N is the total number of observed events, Nj is 
the yield for category j, which corresponds to either 
signal or background, i is the event index and Vj is the 
probability density function (PDF) of the fit component 
j. The PDFs are obtained using MC simulation and are 
modeled as histogram functions. The normalizations 
of single background components (continuum, b — > c, 
and rare) relative to each other arc estimated from the 
MC simulation and fixed in the fit, leading to two free 
parameters in the fit for signal and total background 
yields. Compared to the counting procedure performed 
in the previous analysis [B], this fitting method makes 
use of the signal shape to discriminate between signal 
and background contributions. 

We calculate the significances and the upper limits 
by evaluating the likelihood profile. To take into account 
the systematic uncertainty, we convolve the likelihood 
function with a Gaussian whose width equals the total 
systematic uncertainty. The significance is obtained by 
comparing the likelihood values at maximum and at 



zero signal yield: S 



The upper limit 



on the branching fraction at 90% confidence level (C.L.) 
is evaluated through a Bayesian method by integrating 
the likelihood function from zero to the bound that 
gives 90% of the total area; this assumes a uniform prior 
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FIG. 2: The Eecl distributions for B — > h^vv decays. 
Points with error bars are data; the solid black histogram 
is the total fit result. The blue cross-hatched region is the 
background component; the dashed red histogram shows the 
signal contribution. 



distribution for the branching fraction. We obtain the 
branching fraction using the signal yield iV s i g , the signal 
efficiency e and the total number of BB pairs N B g: 
B = N sig /(e-N BB ). 

To evaluate the sensitivity, simulated experiments 
with the expected amount of background events and 
zero signal events were generated. For each of the 
experiments, an upper limit on the branching fraction 
at 90% C.L. was calculated. The median values of the 
obtained upper limit distributions arc summarized in 
the rightmost column in Table [T] 

The -Eecl distributions in data are shown in Fig. [2] 
superimposed with the fit result. The total numbers of 
observed events, the signal yields, the significances of the 
observed signal, the reconstruction efficiencies and the 
upper limits on the branching fractions are summarized 
in Table [TJ None of the signal modes show a significant 
signal contribution. According to MC studies, the 
enhancements in the K + vD and nvv modes are unlikely 
to be caused by peaking background contributions. The 
signal reconstruction efficiencies are estimated with 
MC simulations using the B — > tv-*' form factors from 
Ref. PR]. 

The systematic uncertainty is dominated by the 
statistical uncertainty of the background model. The 
stringent selection increases the signal to background 
ratio but also reduces the number of MC events in 
the signal box. This leads to a large uncertainty in 
the background shape, despite using an MC sample 
corresponding to five times the data luminosity. To 
estimate the uncertainty, we replace the nominal back- 
ground model with two alternative models compatible 
with the simulation and repeat the fit. The alternative 
background models are Chebyshev polynomials of order 
0, 1 or 2. For each channel, the two models that are most 
compatible with the background distribution are used. 
After the fit with these models, the largest deviation 
of the signal yield from the nominal fit is assigned as 
systematic error, which can vary in size among channels 
due to the different background shapes. The fit bias 
is evaluated through pseudo-experiments with signal 
and background yields set to the observed values. The 
systematic uncertainty due to MC data discrepancy of 
the track and ir° rejection was studied using a 
control sample. Uncertainties associated with the B teig 
reconstruction efficiency, signal MC statistics, particle 
identification, track or particle reconstruction efficiency, 
the total number of the BB pairs and the form factors 
of the signal model are included as well. All systematic 
uncertainties are summarized in Table [Til The total 
systematic uncertainty is calculated by summing all 
contributions in quadrature and is generally smaller 
than the statistical error. 

In conclusion, we have performed a search for 
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TABLE I: Summary of the total number of observed events in the signal box iVtot, the resulting signal yield iV s i g , the significance 
of the observed signal, the signal efficiencies e, the measured and the expected upper limits on the branching fractions at 90% 
C.L. 



Mode 



Significance e, 10 4 Upper limit 



Expected limit 



B+ K + vu 
B° -> K° 3 vv 
B+ -> K*+vv 
B° -> K*°vv 
B + -> n + vv 
B° -> ir°W 
B + p+vv 
B° -> p<W 
B° -> 



43 
4 
21 
10 
107 
6 

90 
31 
3 



+7 



!;e(stat) ± 2.3(syst) 



13.3 

I. 8rJ;|(stat) ± l.O(syst) 
-1.7il;I(stat) ± 1.5(syst) 
-2.3l3° 5 2 (stat) ± 0.9(syst) 

15.2±££(stat) ± 1.4(syst) 
3.5l?;|(stat) ±0.6(syst) 

II. 3±|;|(stat) ± 4.1(syst) 



1.6l 4 . 



"+|;?(stat) ±0.4(syst) 
1.4±jj;jj(stat) ±0.8(syst) 



2. Oct 
0.7a 



2.6a 
1.9a 
1.7a 
OAa 
0.5a 



0.84 
1.47 
1.44 
3.39 
1.66 
1.35 
0.64 
0.58 



< 5.5 x 10~ 5 

< 9.7 x 10~ 5 

< 4.0 x 10~ 5 

< 5.5 x 10~ 5 

< 9.8 x 10~ 5 

< 6.9 x 10~ 5 

< 21.3 x 10 -5 

< 20.8 x 10 -5 

< 12.7 x 10~ 5 



2.2 x 10~ 5 

7.3 x 10" 5 
5.8 x 10" 5 
4.6 x 10~ 5 
3.8 x 10~ 5 

3.6 x 10" 5 
10.2 x 10~ 5 
15.7 x 10~ 5 

8.7 x 10" 5 



TABLE II: Summary of the systematic errors, 
of the signal normalization are given in %. 


The errors 


on the 


signal yield 


are given 


in the number of events and the 


errors 


Channel 


K + vv 


K°vD 


K*+vv 


K*°vD 


TT + VU 


n°vv> 


p + w 




4>vv 


Signal yield [events] 


Background model 


2.1 


0.9 


1.5 


0.5 


0.9 


0.4 


4.0 


0.4 


0.5 


Fit bias 






0.2 


0.6 




0.4 




0.1 


0.6 


Signal normalization /%/ 


Track and tv° rejection 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


Stag correction 


4.2 


4.5 


4.2 


4.5 


4.2 


4.5 


4.2 


4.5 


4.5 


Signal MC statistics 


1.2 


3.5 


3.7 


2.8 


1.5 


2.1 


2.3 


3.3 


2.6 


Track, 7r° and Kg reconstruction efficiency 


0.3 


2.3 


4.1 


0.4 


0.4 


4.0 


4.2 


0.7 


1.4 


Particle identification 


2.0 


4.0 


2.0 


4.0 


2.0 




2.0 


4.0 


4.0 


N B b 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


Form factors 


2.0 


5.4 


3.8 


6.4 


1.9 


1.6 


2.9 


4.5 


7.5 



B — > hy>vi> decays in nine different modes with a fully 
reconstructed B tag on a data sample of 772 x 10 6 BB 
pairs collected at the T(4S) resonance with the Belle 
detector. No significant signal is observed and we set 
upper limits on the branching fraction at 90% C.L. The 
results of this analysis supersede the previous results 
from Belle 6 . The limits reported here for K* + vv, 
■n^vv, tt {) vv and p°vv are the most stringent constraints 
to date [16]. These limits are above SM predictions and 
allow room for new physics contributions. The upcom- 
ing Belle II experiment [T7] should be able to reach a 
sensitivity high enough to probe the SM predictions for 
exclusive b — > svv decays. 
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